

O 



INVESTOR IN PEOPLE 



The Patent Office 



Concept House 
Cardiff Road 
Newport — . 




I, the undersigned, being an officer duly authorised in accordance with Section 74(1) and (4) 
of the Deregulation & Contracting Out Act 1994, to sign and issue certificates on behalf of the 
Comptroller-General, hereby certify that annexed hereto is a true copy of the documents as 
originally filed in connection with the patent application identified therein. 



In accordance with the Patents (Companies Re-registration) Rules 1982, if a company named 
in this certificate and any accompanying documents has re-registered under the Companies Act 
1980 with the same name as that with which it was registered immediately before re- 
registration save for the substitution as, or inclusion as, the last part of the name of the words 
"public limited company" or their equivalents in Welsh, references to the name of the company 
in this certificate and any accompanying documents shall be treated as references to the name 
with which it is so re-registered. 

In accordance with the rules, the words "public limited company" may be replaced by p. I.e., 
pic, P.L.C. or PLC. 



Re-registration under the Companies Act does not constitute a new legal entity but merely 
subjects the company to certain additional company law rules. 



Signed VA KjZ^M,i 

Dated 30 Ochobtr dlOOO 




PRIORITY 
DOCUMENT 



SUBMITTED OR TRANSMITTED IN 
COMPLIANCE WITH RULE 17.Ua) OR (b) 



An Executive Agency of the Department of Trade and Industry 



itents Form 1/77 



Qfl5oe 




Request for grant of a patent 

(See the notes on the back of ibis form. You con aiso get 
an explanatory leaflet from the Patent Office to bttip 
youfM in this form} 




The Patent OfQce 

Cardiff Road 
Newport 
Gwent NP9 IRH 



1 . Your reference 



140058 



210CT99 E^85665~l D02ilL 



P01/77QC 0.00 - 992^g:3.S 



2. Patent application number / /? o -> 

ObePatem Office turn fttt in tMs part) QQ n R A A/- 6 O 3 " 5^ 



3. Full name, address and postcode of the or of 

each applicant (underline aU surnames) 

Patents ADP number df^ou knoiv U) l)L 

If the applicant is a coiporate body, give me 
country/state of its incorporation 



BTG INTERNATIONAL LIMITED 

10 Fleet Place 

Limeburner Lane 

London 

EC4M 7SB 

England 

G6 



4. Utie of the invention 



PERMANENT MAGNET AND SHIM DESIGN 



5. Name of your agent Of you have one) 

"Address for service" in the United Kingdom 
to which all correspondence should be sent 
(induding the postcode) 



Patents ADP number (if you knotv it) 



SCOTT, Susan Marg^et 

BTG INTERNATIONAL LIMITED 

10 Fleet Place 

Limeburner Lane 

London 

EC4M 7SB 

England 



6. If you arc declaring priority from one or more 
earUcr patent applications, give the country 
' and the date of filing of the or of each of these 
earlier applications and (if you know it) the or 

each application number 



Country Priority application iiunU>er 

(if you know it) 



Date of filing 
fdkxy / month /year) 



7, If this application is divided or otherwise 
derived from an earlier UK application, 
give the niimber and the filing date of 
the earlier application 



Number of earticr :q;>plication 



Date of filing 
(day / month /year) 



Is a Statement of inventorship and of right 
to grant of a patent required in support of 
this request? (Answer -Yes' if: 

a) any applicant named in part 3 is not an inventor, or 

b) Ufere is an intfentor who is ma named as an 
€ippUcant, or 

c) any named applicant is a corpontte body. 
See note (d)) 



YES 



Patents Form 1/77 



u 



Patents Form 1/77 



9. 



Enter the number of slU^for any of the 
following items you are filing with this form. 
Do not count copies of the same document 

Continuation sheets of this form 
Description 

Abstract 




18 

3 
1 




10. If you arc also filing any of the following, 
state how many against each item. 

Priority documents 
Translations of priority documents 



Statement of inventorship and right 

to grant of a patent {Patents Form 7/77) 



Request for preliminary examination 

and search (patents Form 9/77) 

Request for substantive examination 
(Patents Form 10/77) 



Any other documents 

Cpleose ^>ecify) 



1 1 . I/We request the grant of a patent on the basis of this application. 

Signature ^ ^^"^ . Date 20,10.99 



SCOTT, Susan Margret 20 October 1999 



12. Name and daytime telephone number of - ii«Tr*.n-« .«« --i- e-iw^ if?«ji 

pc«on to coAact in the united Kingdom "^^^CHER 171 575-1584 ^ 

Warning 

AJter an application for a patent has been filed, the Comptroller of the Patent Office will consider whether publication 
or communication of the invention should be prohibited or restricted under Section 22 of the Patents Act 1977. You 
Witt be informed if it is necessary to prohibit or restrict your invention in this way. Furthermore, if you live in the 
United Kingdom, Section 23 of the Patents Act 1977 stops you from applying for a patent abroad without first getting 
written permission from the Patent Office unless an application has been filed at least 6 weeks beforehand in the 
United Kingdom for a patent for the same invention and either no direction prohibiting publication or 
communication has been given, or any such direction has been revoked 



Notes 

a) If you need help tofiU in this form or you have any questions, please contact the Patent Office on 064$ 500505. 
h) Write your answers in capital letters using black ink or you may type them. 

c) If there is not enough space for all the relevant details on any part of this form, please continue on a separate 
sheet of paper and write "see continuation sheet" in the relevant part(sX Any continuaHon sheet should be 
attached to this form. 

d) If you have answered 'Yes' Patents Form 7/77 will need to be filed 

e) Once you have filled in the form you must remember to sign and date it. 

f) For details of the fee and ways to pay please contact the Patent Office. 



P^enls Form 1/77 



. F:\140\058\GB01SPEaASFILEDJXX; 

r 

PERMANENT MAGNET AND SHIM DESIGN 

FIELD OF THE INVENTION 
This invention relates to the design of permanent magnetic structures for 
5 generating a specific magnetic field in a Region Of Interest (ROI) or modifying an 
existing magnetic field in an ROI using shims. 

BACKGROUND OF THE INVENTION 
Permanent magnets are of great value in many fields of science and 
engineering. In particular, they find application in the field of Magnetic Resonance 
10 Imaging (MRI) and Spectroscopy (MRS). Within these latter applications it is of 

paramount importance that extremely high degrees of magnetic field homogeneity 
(of the order of one part per million) are achieved within an ROI. 
In permanent magnets the distribution of the magnetic material determines 
how the resulting magnetic field varies over space. In both permanent magnetic and 
15 electromagnetic MRI and MRS machines, manufacturing imperfections result in 
unacceptable inhomogeneities. To improve the homogeneity of an existing field 
over an ROI it is ioiown to place permanent magnetic materials (termed passive 
shims) in the magnetic field to achieve this improvement of homogeneity. 

The problem of distributing magnetic material to produce, or modify an 
20 existing magnetic field is one . of a class known collectively as the near field 
problem. In particular, it requires the determination of the distribution of material 
to produce a field of particular character - varying in a specific manner in space, for 
example a uniform field in a proximate region. The total magnetic field from a 
magnetic structure is a nonlinear convolution of the spatial distribution of material 
25 with the field of a magnetic dipole. While it is relatively straightforward to 
determine the resulting field from a given distribution of magnetic material, 
generally, the components of a convolution are not invertable. Consequently, 
detemiining the optimum distribution of magnetic material to produce a magnetic 
field of a specific type in an ROI has proven to be extremely difficult. 
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One known method of electromagnetic coil design is the Target field approach 
as disclosed in European Patent No. EP-A- 0 252 634 or its US equivalent No. 
4896129. As it is based on an assumed a priori mathematical expression for the 
total field, an optimum resulting field is hard to achieve. 
5 An alternative approach is disclosed in European Patent EP-A-0601 101 and its 

US equivalent No. 5266913, the content of which is hereby incorporated by 
reference. It is concerned with the design of shielded electromagnetic coils, and 
employs forward substitution and least squares minimization to optimize the field 
produced by electromagnets specifically within the ROI, whilst achieving a null 
10 field outside of the confines of the coil structure. 

There are significant differences in the form of the magnetic fields generated 
by active electromagnets, and those resulting fi*om the use of magnetic materials. 
Consequently, electromagnetic coil design and the design of permanent magnets 
involve quite separate design considerations. 
15 The object of the present invention is to achieve an improved method of 

designing permanent magnetic structures. 

SUMMARY OF THE INVENTION 
In accordance with the present invention there is provided a method of 
designing a permanent magnetic structure for generating a desired magnetic field in a 
20 region of interest (ROI) comprising: 

-setting the dimensions of the permanent magnetic structure 
-setting a hypothetical distribution of magnetic material over the surface of 
said structure for each member of a finite set of spatial frequency harmonics of unit 
magnetisation amplitude 
25 -calculating the magnetic field in the ROI from each member of said set 

-calculating amplitudes of said harmonics by the method of least squares in 
order to achieve the optimum s^proximation to the desired field in the ROI 

-scaling the members of said set according to said optimised harmonic 
amplitudes 
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-establishing, by summing said scaled members, the required distribution of 
magnetic material in order to generate the desired magnetic field. 

The permanent me^netic structure may either be a permanent ms^et 
assembly for producing a magnetic field ab initio, or may be a permanent magnetic 
5 structure for modifying an existing magnetic field, by means of shimming. 

The method may also be applied iteratively in order to take into account the 
fuiite thickness of the required distribution magnetic material, which will inevitably 
have influence on the field. 

Further preferred embodiments of the present invention are presented in the 
10 dependent claims. 

BRIEF DESCRIPTION OF THE DIAGRAMS 
Figure 1 shows a schematic representation two parallel, finite, flat plates upon 
which magnetic material is distributed. 

Figure 2 shows an alternative distribution composed of a hollow elliptical, 
1 5 cylindrical structure upon which magnetic material is distributed. 

Figure 3 shows the magnetic dipolar field (Bd) produced by magnetic material, 
resolved into orthogonal components Bdr , Bdz and Bda. 

Figure 4 shows schematically a hollow, cylindrical structure with a circular cross- 
section, upon which magnetic material is distributed. 
20 Figure 5 shows the permissible, even, cosinusoidal axial harmonics of the magnetic 
material distribution on a structure of finite length. 

Figure 6 is a plot of the axial mc^etic field Brz i^ainst distance z resultix^ firom an 
infinitesimal circular ring of uniformly magnetic material of unit magnetization 
calculated directly by numerical integration. 
25 Figure 7 shows a plot of the Fourier Transform (FT) in the axial direction (z) of the 
axial magnetic field of Figure 6. 

Figure 8 shows a plot of the FT, in the axial direction, of the axial magnetic field 
resulting from an infinitesimal circular ring of uniformly magnetic material of unit 
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magnetization evaluated by an alternative, analytical method, directly in the Fourier 
Kernel domain (k-space). 

Figure 9 shows a plot of the axial distribution of magnetic material that >vill achieve 
an optimal, uniform field, designed in accordance with a first embodiment of the 
5 present invention. 

Figure 10 shows a schematic representation of a magnetic structure designed in 
accordance with the invention, to generate the distribution shown in Figure 9. 
Figure 1 1 shows a plot of the resulting axial magnetic field in the axial direction 
produced by the distribution of magnetic material shown in Figure 9. 
10 Figure 12 shows a plot of the deviation of the above magnetic field from the 
desired, uniform one, in parts per million (ppm), over the ROI. 
Figure 13 shows a plot of the deviation of the above field over, and proximate to, 
the ROI. 

Figure 14 shows how a distribution of paramagnetic and diamagnetic materials on a 
15 surface can be azranged to produce a sinusoidal variation of magnetization along 
that surface. 

Figure 15 shows a plot of the amplitude of the optimum axial distribution of 
magnetic material over a cylindrical structure to provide a magnetic field that varies 
linearly in a transverse direction in the ROI, designed in accordance with the present 
20 invention. 

Figure 16 shows a schematic representation of how a shim magnetic structure 
designed in accordance with the invention and using diamagnetic and paramagnetic 
materials as illustrated in Figure 14 could be employed to generate the distribution 
shown in Figure 15. 

25 Figure 17 shows a plot of the resulting axial component of the gradient magnetic 
field at a series of radii in the axial direction produced by the distribution of 
magnetic material shown in Figure 15. 
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Figure 1 8 shows how the sinusoidal variation of thickness in a single unidirectional 
magnetic material can provide a sinusoidal variation of magnetization for the 
purpose of shimming. 
Theory. 

5 The design of ms^ets for MRI and MRS is commonly based upon the forms 

of simple geometric surfaces. The two most common of these are: a structure 
composed of parallel, finite, flat surfaces between which a uniform magnetic field is 
created as illustrated in Figure 1; and a structure with an elliptical, cylindrical 
opening within which a uniform field is produced by the surrounding distribution of 
1 0 material as illustrated in Figure 2, 

With reference to Figure 3, the magnetic dipolar vector field (Bd) produced by 
magnetic material can be resolved into its cylindrical polar components Bdz, Bde and 
Bdr. In these coordinates (r, 0, z), the vector field of such an infinitesimal dipole 
has the general form: 

15 

Bdz{r,z) - ^ m (22? . i^) / 47C (r^ + 2?)^'^ (la) 
Bde(r.z)-0 (lb) 

and 

20 Bd^r,z) = n m 3zr / 4it (r^ + 2?)^^ (Ic) 

where \i is the permeability of the medium in the ROI and m is the magnetic dipole 
moment. The field due to the circulation of an electron in an atom approximates 
this very closely for distances (r^+ z^)''^ that are very much greater that the radius 
25 of the orbit (typically 10"'*^ m). For the purposes of passive magnet design within 
the present context this approximation holds extremely well. 

The magnetic field from matter in bulk can be derived firom equation 1. For 
the axial component of the magnetic field this is given by: 
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B,(r,z) = f {Mz(x,y,z) (2(2-2')^ - [(x-xf + (y-yf ]) / 47t|r - ff} dx' (2) 

where the underlining of individual items indicates that they are vector quantities 
and Mz is the axial component of magnetization (M) . Similar expressions can be 
5 derived for the field in the other directions. 

For paramagnetic and diamagnetic materials the magnetization M is given by: 

M = X„.Bm (3) 

10 where Bm is the main, magnetizing field and Xm is the susceptibility. In general, 
Xm is a tensor of rank two. However, for isotropic materials it reduces to one of 
Tank zero, e.g., for diamagnetic copper Xm = -l.O x 10*^; and for paramagnetic 
aluminium Xm = 2.3 x 10'^. 

Equation (2) shows that the magnetic field is the convolution of the 

15 magnetization distribution with a spadal dispersion fimction. Problems of this 
nature can be conveniently solved by Fourier transforming into the kernel space (k- 
space) where the convolution becomes a simple multiplication, which greatly 
facilitates subsequent calculations. Indeed, for certain simple geometric structures, 
such as cylinders, there is no simple analytic expression for the magnetic field in 

20 real space, while there is one for the field in k-space. The existence of such an 
expression in k-space expedites the process of handling the above equations but it is 
not a necessaiy condition of the design method of the present invention. Indeed, the 
method is generally applicable, as is illustrated in the following examples. 

25 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The design of a permanent magnetic structure according to the invention will 
now be described in general terms. The geometry and overall dimensions of the 
desired structure are set, as is a hypothetical distribution of magnetic material over 
the surface of the structure for each member of a finite set of spatial firequency 
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harmonics of unit magnetization amplitude. The field within the ROI fix)m an 
element of the ms^netic material of unit magnetization can then be determined. 
When numerical methods are employed the field must be determined to a point 
where it £q}proximates zero in the directions in which it is to be Fourier Transformed 
5 (FTed). This field is then FTed. The resulting field in Fourier kernel space (k- 
space) is multiplied by the FT's of a set of the permissible harmonics (the lowest 
spatial frequency harmonics of the initial magnetic field distribution) of the 
magnetization on the surface. These resulting fields in k-space are subsequently 
inverse FPed to yield the fields in real space that are produced by each of the 
10 permissible harmonics in the selected set. These fields are fitted to the required 
field within the ROI to determine the amplitudes of the magnetization harmonics 
that minimize the sum of the squares of the deviation of the actual field fi-om the 
required one. Fitting is performed using sets containing successively higher order 
harmonic components until the sum of the squares of the deviations of the 
15 achievable field firom the desired one becomes sufficiently small. Eventually a 
point will be reached where deviations due to manufacturing tolerances will 
dominate, and no further practical improvement will be possible. 

For simple geometric structures, such as a circular cylinder, this process can 
be facilitated by directly calculating the fields in Fourier k-space. 
20 The advantage of least squares fitting lies in its convenient mathematical 

properties. For example, it commutes with the Fourier transform; as the FT itself is 
a specific example of least squares minimization. 

The following discussion serves to illustrate the method of the present 
invention only and is not intended to limit its general applicability. 
25 In a first embodiment of the present invention, the method is applied to design 

a magnet where magnetic material is distributed on the curved surface of a cylinder 
to produce a uniform field. 

For the present purposes it will be supposed that it is required to create a 
uniform axial magnetic field Bz of 0.1 T over a cylindrical ROI of length 0.4m and 
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radius 0.125m by distributing magnetic material on the curved surface of a circular 
cylinder, such as that shown in Figure 4, having a length L of 2m and radius a of 
0.25m. The cylindrical ROI and surrounding cylindrical structure are concentric 
and coaxial along the z axis. For simplicity it will be assumed that the thickness of 
5 the magnetic shell is small compared to the diameter of the magnet, for the present 
purposes infinitesimal. In practice, a more elaborate calculation may be required to 
take into account the thickness, as described later. 

Conveniently, optimization is performed over a representative set of radii 
spanning the ROI, e.g. six radii at intervals of 20% from 0 to 100% of the radius of 

10 the ROI within the ROI. However, for the present example, optimization over just 
two curved cylindrical surfaces of diameter of Om and 0.25 m and lengths 0.4m, was 
sufficient, as demonstrated by table 3 and Figure 12 (explained below). 

The distribution of magnetic material on the sur&ce of the structure shown 
in figure 4 can be conveniendy described in terms of a series of orthogonal 

15 azimuthal and axial, sinusoidal and cosinusoidal, harmonics. The azimuthal 
harmonics form a disciete set of cosinusoidal terms cos (m<|> + a) and the axial 
harmonics form a truncated set of H(z) cos (n9 + P), as illustrated in Figure 5, 
where H(z) is the unit step function given by: 

20 for|z| ^L/2 H(z)=l 

(4) 

for|z| ^L/2 H(z) = 0 

where L is the length of the cylindrical surface and a and P are phase terms that take 
25 account of both the sinusoidal and cosinusoidal terms. 

In many applications the set of harmonics is severely curtailed by 

considerations of synunetry. For example, for the present uniform field there is no 

azimuthal variation and all azimuthal terms other than m=0 are zero. 

Consequently, the problem of optimization reduces to a one-dimensional one of 
30 determining just the axied variation. Further, as the ROI and the structure are 



-8- 



F:M40V058VGB0I SPECVASFILED.DOC 



concentric the resulting even symmetry of the system implies that all sinusoidal 
terms are also zero leaving only the purely cosinusoidal terms. It remains to 
determine the amplitudes of each of the significant cosinusoidal terms to produce an 
optimum uniformity of the field within the ROI on the basis of the minimum sum of 
least squares deviations of the field. 

Now the field from an elemental ms^etic dipolar ring of unit magnetization 
can be determined by numerically integrating equation (2). Figure 6 shows how the 
resulting magnetic field varies with distance, in the axial direction, at a radius of 0.1 
m from such a dipolar ring of diameter 0.5 m. 

From equation (2) we also know that the magnetic field fix>m individual axial 
harmonics is the convolution of the field from this ring with the harmonic. This 
convolution is most conveniently performed by simple multiplication in the Fourier 
domain. Consequently, the next step in the design process is to FT this field in the 
axial direction z; the resulting FT is shown in Figure 7. 

Alternatively, for a simple geometric surface such as this, the FT in the axial 
direction of the axial component Brz of the ring field, Br, can be expressed 
analytically as: 



where Im and Km are the modified Bessel functions. 

Figure 8 shows the FT in the axial direction of the axial component Brz of the 
ring field, Br, evaluated directly by equation 5a. This is substantially identical to the 
field derived by numerically integrating the generally applicable equation (2) and 
then Fourier Transforming it, as can be seen by comparing Figures 7 and 8. 

Although evaluating the FT directly by means of equation (5a) facilitates the 
process, this method is not applicable to structures where an exact solution in the 
Fourier kernel domain is not available. The longer route of numerical integrating 



f or r > a 



for r ^ a 



Brz (r.m,k) = - pi a U|k| r) Km(|k|a) / 2n 
Brz (r,m,k) - - ^ a I„(lk| a) Kn,(lk|r) / In 



(5a) 
(5b) 
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equation (2) and subsequent FT'ing has general application for any structure, as it 
does not rely upon there being an analytic form of the equations in either real or kr 
space. 

Now, for an open cylinder of restricted length, L, the total magnetization, Mz, 
can be expressed as the sum of the pennissible axial harmonic modes, n, on the 
surface as follows: 

Mz(a, z) = In C„ H(z) cos (2 71 n z / L - <t» ) (6) 

where Cn is the amplitude of the nth magnetization harmonic. The FT of this in the 
z-direction is: 

MjHsl, k) = ir L Si, C„ [e^ sinc(k L/2 - nn) + e*** sinc(k L/2 + mc)] (7) 



Therefore, the FT of the total resultant field &om the cylindrical sur&ce is 
simply the product of the FT of the magnetization with the FT of the elementary 
ring field within this domain, i.e.: 

Bz(r, m, k) = Mz(a, k) Br^ (r, m, k) (8) 

Optimization of the distribution can is then solved by least squares 
minimization of the first few terms (Cn) over the ROI in real, as against kernel, 
space. The number of terms is dependent upon the geometries of the coil structure 
and the ROI. This can be determined by solving with increasing numbers of terms 
(Cn) until a sufficiently homogeneous field has been attained. 

In the above example the first 3, 4 and 5 harmonic amplitudes (coefficients) 
were employed and the worst deviation from the field in the ROI determined (see 
Table 1). Optimization over more than 5 harmonic coefficients has no further 
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significant effect. As can be seen from Table 1, the maximum deviation of the field 
in the ROI drops to ~ Ippm within the first 5 terms. In any case, practical 
considerations, such as manufacturing tolerances and the finite thickness of the 
structure, mean that further improvements in field uniformity are impractical. 
5 Further approximation can be achieved by shimming the magnet in some fashion, 
for example along the lines discussed in the second embodiment below. 



Table 1: Field deviation with number of harmonic components employed. 



No. of Harmonics 


Range of Harmonics 


Max deviation in ROI 


3 


0-2 


200 ppm 


4 


0-3 


16 ppm 


5 


0-4 


1 ppm 



10 The harmonic amplimdes of table 2 were used to determine the axial 

distribution of magnetic material (M = 1.28 x 10* amp/m, equivalent to the 
magnetization of saturated iron) on the surface of the hollow cylindrical structure, as 
shown in Figure 9. 

Figure 10 shows a practical implementation of a magnetic structure designed 
1 5 in accordance with the invention, to generate the distribution shown in Figure 9. 



Table 2: The harmonic amplitudes for a uniform field. 



Harmonic No. 


Amplitude (amps/m) 


0 


-1,143,966 


1 


-220,943 


2 


37,362 


3 


-6,821 


4 


749 
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Figure 1 1 shows the resulting magnetic field at a series of radii (0, 2.5, 5, 7.5, 
10 and 12.5 cm) plotted in the axial direction produced by the distribution of 
magnetic material, according to the present invention, as shown in Figure 9. Figure 
12 shows the deviation of the above magnetic field fix>m the desired, uniform one, 
in parts per million (ppm), over the ROI. Figure 13 shows how the field diveiges 
rapidly outside die ROI in the axial direction, but remains reasonably homogeneous 
in the radial direction; see also Table 3 below. 



Table 3: Maximiun field deviation in the ROI at each radius. 



Radius cm (% fraction of magnet) 


Maximum deviation in ROI 


0.0 (0) 


0.460 ppm 


2.5 (10) 


0.160 ppm 


5.0 (20) 


-0.541 ppm 


7.5 (30) 


-1.133 ppm 


10.0 (40) 


-1.011 ppm 


12.5 (50) 


0.164 ppm 


15.0 (60) 


2.01 1 ppm 


20.0 (70) 


-4.730 ppm 



As demonstrated by the above figiires, the resulting field is uniform to a very high 
degree in the central region. 

Figtire 1 1 serves to highlight the fimdamentai inapplicability of methods that 
attempt to invert the field equations to the design of permanent ms^etic structures. 
Due to the complexity of the field, it is inconceivable that even someone highly 
conversant with the design of magnets could arbitrarily specify a field such as that 
of Figure 1 1 with the required degree of accuracy to ensure a practically realizable 
distribution that would yield a reasonably good approximation to the desired field in 
the ROI; essentially this would require a knowledge of the result of the process 
before the process had begun. Certainly, this would not yield the uniformity attained 
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by the present technique which only requires the specification of the known, desired 
field within the ROI. 

In accordance with a second embodiment of the present invention, the design 
method is applied to homogenize an existing field by distributing magnetic material 
on the curved surface of a circular cylinder, of the type shown in Figure 4. For the 
purpose of illustrating this example, the existing field is assumed to be oriented in 
the z-direction and is substantially uniform except for a first order (linear) 
transverse variation (inhombgeneity) in the ROI. It is to be appreciated that the 
design method is by no means limited to compensating for a linear variation, but is 
generally applicable to shimming for any inhomogeneity. 

As in the previous embodiment, the magnetic surface is arbitrarily assumed to 
be 0.5 m in diameter and 2 m in length, with a cylindrical ROI of diameter 0.25 m, 
and length 0.4 m. Optimization was performed on a representative curved 
cylindrical surface of diameter 0.2 m and length 0.4m to produce a transverse linear 
gradient field of 0.1 milliTesla per metre. Again, table 5 below and Figure 17 
demonstrate that this is indeed representative of the ROI. 

For simplicity it will again be assumed that the thickness of the magnetic shell 
is small compared to the diameter of the magnet, for the present puiposes 
infinitesimal. In practice, for the creation of shim structures this approximation will 
hold and more elaborate calculations will generally not be necessary. 

From symmetry considerations, an azimuthal variation of cos (<|i), 
corresponding to a linear transverse variation, in the field can only be produced by 
a first order variation (m=l) in the magnetic material. This again reduces the 
problem to a one-dimensional one of determining just the axial variation. 

The FT of the field can be determined directly from equation 5 a. Optimization 
of the shim is then conveniently solved by least squares minimization of the first 
few terms (Cn) over the ROI in real, as against kernel, space. 

In this embodiment coefficients 1, 2, 3 and 4 were necessary as can be seen 
from table 4. From inspection, it was noted that the zeroth coefficient contributed 
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negligibly to the central field, so this too could be omitted. The amplitude of the 
required axial distribution is shown in Figure 15. The distribution varies azimuthally 
as cos(0) (with amplitude = 1.28 x 10^ amp/m) on the surface of the cylindrical 
structure to provide a magnetic field that varies optimally, linearly, transversely 
(Bz(x) = X.100 microTesla/m) over the ROL 

Figure 16 shows how a shim magnetic structure designed in accordance with 
the invention and using diamagnetic and paramagnetic materials as illustrated in 
Figure 14 may be implemented in practice, to generate the distribution shown in 
Figure 15. 

Figure 17 shows the resulting gradient magnetic field at a series of radii (0, 
2.5, 5, 7.5, 10 and 12.5 cm, 15 cm) , plotted in the axial direction, produced by the 
distribution of magnetic material shown in Figure 15. Note that the linearity remains 
excellent far outside the ROL 



Table 4: The harmonic amplitudes for a linear transverse field. 



Harmonic (fo. 


Amplitude (amps/m) 


1 


247 


2 


-59 


3 


15 


4 


-2 



The resulting magnetic field is linear to a very high degree in the central 
region, see table 5 and figure 17. Certamly, the degree of linearity is well within the 
limitations of first order approximation adopted and of the errors associated with 
manufacturing. Consequently, such a shim designed in accordance with the present 
invention will correct the linear inhomogeneity to a very high order of accuracy. As 
mentioned above, the design method of the present invention is not limited to 
correction of linear variations, and a similar to degree of correction can be 
anticipated generally, for any conceivable inhomogeneity. 
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Table 5: Center field in the RQI at each radius for a linear transverse shim. 



Radius cm (% fraction of magnet) 


Bz (microTesIa) 


0.0 (0) 


0.000000 


2.5 (10) 


2.499998 


5.0 (20) 


4.999995 


7.5 (30) 


7.499989 


10.0 (40) 


9.999974 


12.5 (50) 


12.499929 



It is to be appreciated that, generally, the lowest spatial frequency harmonics 
are chosen, but this will not always be necessary as witnessed by the example 
above. 

The number of terms over which the optimization is performed may be 
reduced by incorporation of an additional, hypothetical, arbitrary distribution of 
magnetic material over the surface of the structure into the set of spatial frequency 
harmonics. By ''arbitrary'* it is to be understood that a suitable distribution is 
determined by appropriate consideration of the desired field. For example, to 
achieve a field that varies linearly along a flat sur&ce of magnetic material, an 
appropriate distribution might be one whose magnetization also varies lineariy 
along that surface. 

The method of the invention can be applied to the design of shims in several 
ways. One of these is to fit the dipole harmonics directly to the entire field 
inhomogeneities, as illustrated above. Another is to fit the dipole harmonics to 
individual spherical harmonic terms of the field. The rationale for the latter is that it 
is common practice to map the field in these linearly independent spherical 
harmonic terms. The required shim is then the linear combinations of the shims for 
the individual spherical harmonic terms weighted appropriately. This second 
method is advantageous when shimming a series of magnets of the same physical 
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dimensions, as the shim dimensions remain constant and consequently the spherical 
harmonic shims form a standard set. The resulting shim from the example above 
effectively also represents one component of such a set 

The required distribution of magnetization can be achieved in several ways. In 
one preferred embodiment the thickness of a single magnetic material can be varied 
in proportion to the required distribution. This is perfectly adequate for the 
production of uniform fields such as described in the first embodiment provided the 
thickness of the material is small in comparison with the sizes of the structure and 
the ROI. Where this is not the case the process must be refined by shimming along 
the lines outlined for the second embodiment 

For the purpose of shinuning the use of a single imidirectional material cannot 
produce a magnetization that varies sinusoidally, but if an additional uniform field 
can be tolerated simply varying the thickness of the material sinusoidally in the 
manner shown in Figure 18 produces acceptable permissible harmonic 
magnetization components of the form A-cos(m0) and A-sin(m9), in the example 
illustrated, 2-sin(39). 

In another preferred embodiment distributing materials of dififering 
magnetization in the vicinity of the ROI can be used to achieve the desired local 
magnetization; By using differing thicknesses of materials that have opposite 
susceptibilities harmonic components of the form cos(m6) and sin(m6) are possible. 
For example, for a cylindrical surface consisting of a semicircular (0 < 6 < tc) 
distribution of paramagnetic aluminium (Xm - 2.3x10'^), whose thickness varies as 
0.435 A |sin{0)| and an opposing semicircular (-Jt<9<0) distribution of 
diamagnetic copper (xm = -10x10'^) whose thickness varies as A |sin(0)|, produces a 
magnetization that varies as A sin(9) to first order over the whole surface 
(-71 <9 <7r), in a similar manner to that shown in Figiu^ 14. Alternatively, the 
magnetization distribution shown in Figure 15 could be achieved by machining the 
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above sinusoidal contour on cylinders of paramagnetic and diamagnetic material in 
this manner. 

A material of any required magnetism can be constructed by combining 
paramagnetic and diamagnetic materials in the required proportions. For example, a 
5 material that has effectively zero magnetization can be formed &om a combination 
of 30.3% aluminium and 69.7% copper. This could be a great value in 
manufacturing objects for use in magnetic fields that would not affect the field, or 
be affected by the field. In practice, the proportions may have to be altered slighdy 
to take into account the higher order effects when such materials are combined. 
10 For either creating a uniform field or altering an existing one, the 

thickness of the magnet will be a significant firaction of the diameter and the 
consequent placement of material will affect the actual attainable field. To take 
account of this, the method can be applied iteratively as follows, by: 

-calculating the distribution of magnetic material to produce an optimum 
15 approximation to a desired magnetic field in the ROI, in accordance with the 
method described herein 

-detemiining the total resultant field in the ROI taking into account the 
thicknesses and placements of magnetic materials involved 

-calculating a corrected magnetic material distribution representative of the 
20 difference between said total resultant field and said desired field 

-repeating the latter two steps until the field attained approximates the 
desired one in the ROI to a sufficient degree. 

It is to be appreciated that the homogeneity of the field generated by a 
magnetic structure designed in accordance with the present invention can be further 
25 enhanced by additional active shimming. Resistive electromagnetic shimming may 
subsequently be employed to further improve the homogeneity during the use of the 
magnet. This is particularly relevant at high field strengths with gradient echo 
techniques where the natural susceptibility of the sample distorts the field, in other 
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# 

words, where resistive electromagnetic shimming can be employed to improve the 
homogeneity for each sample. 

Whilst the above embodiments have been described in the context of magnets 
of cylindrical form, it is appreciated that the method of the invention is equally 
applicable to the design of magnets of other geometric forms, one example of which 
is the parallel flat plate arrangement shown in Figure 1. 



- 18- 



F:\ 1 40\058\GBO I SPEC\ASFILED.DOC 




CLAIMS. 

1 . A method of designing a permanent magnetic structure for generating a 

desired magnetic field in a region of interest (ROI) comprising: 

-setting the dimensions of the permanent magnetic structure 

-setting a hypothetical distribution of magnetic material over the surface of 

said structure for each member of a finite set of spatial firequency harmonics of unit 

magnetisation amplitude 

-calculating the magnetic field in the ROI fi-om each member of said set 
-calculating amplitudes of said harmonics by the method of least squares in 

order to achieve the optimum approximation to the desired field in the ROI 

-scaling the members of said set according to said optimised harmonic 

amplitudes 

-establishing, by summing said scaled members, the required distribution of 
magnetic material in order to generate the desired magnetic field. 

2. A method according to claim 1 wherein the desired field represents a 
modification of an existing magnetic field. 

3. A method according to claim 2 wherein the desured field represents the 
deviation fi-om a required field generated by an existing magnetic structure. 

4. A method according to claim 3 wherein the desired field represents one of 
the harmonics of the magnetic field over the ROI. 

5. A method according to claim 4 wherein each said harmonic results in a 
shim, to create a set of shims representative of the weighted sum of the individual 
harmonics. 

6. A method according to any of claims 1 to 5 wherein the desired field is 
produced by a distribution of magnetic material over one or more flat surfaces. 

7. A method according to any of claims 1 to 5 wherein the desired field is 
produced by a distribution of magnetic material over an elliptical cylindrical shell. 
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8. A method according to any of claims 1 to 5 wherein the desired field is 
produced by a distribution of magnetic material over the surface of a circular 
cylindrical shell. 

9. A method according to claim 7 or claim 8 wherein the thickness of 
5 magnetic material is varied in an azimuthal direction in accordance with the desired 

field. 

10. A method according to claim 1 wherein the desired field is produced by a 
combination of paramagnetic and diamagnetic materials. 

11. A method of designing a permanent magnetic structure for generating a 
1 0 desired magnetic field in a region of interest (ROI) on the basis of the minimization of 

the sum of the least squares of the deviations over the ROI. 

12. A method according to any of claims 1-11 wherein an additional, 
hypothetical, arbitrary distribution of magnetic material over said surface is 
incoiporated into said set 

15 13. A method according to claims 2 or 3 comprising: 

-calculating the distribution of magnetic material to produce an optimum 
approximation to a desired magnetic field in the ROI, in accordance with the method 
described herein 

-determining the total resultant field in the ROI taking into account the 
20 thicknesses and placements of magnetic materials involved 

-calculating a coirected magnetic material distribution representative of the 
difference between said total resultant field and said desired field 

-repeating the latter two steps until the field attained approximates the desired 
one in the ROI to a sufficient degree. 
25 14. A method according to any of claims 1 to 5 where the magnetic fields of 

said set of harmonics are evaluated by an actual step involving numerical integration 
in real space. 



-20- 



F:\I40\OS8\CBOISPEC\ASFILED.DOC 



m ■ ■ 

15. A method according to any of claims 1 to 5 where the Fourier 
transformers of the magnetic fields of said set of harmonics are evaluated directly in 
Fourier Kernel space. 

16. A permanent magnetic structure designed in accordance with the method 
of any of claims 1 to 15. 
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ABSTRACT 

PERMANENT MAGNET AND SHIM DESIGN 

A general method for the design of structures composed of permanent 
magnetic material for producing and modifying magnetic fields. The method 
employs a finite set of permissible spatial harmonic modes upon a surface. In 
particular, it has applications were the structures are finite and open, for which an 
exact solution is not possible. The method determines an optimum distribution 
based upon a least squares minimization of the difference between the achievable 
and desired field within a Region Of Interest (ROI). The method also has 
application in passively improving the homogeneity of existing magnetic fields 
(shimming) by distributing magnetic materials in the vicinity of the ROI. 
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FT of Bz versus Distance. 
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